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Comparative analysis of different types of initial mass spectra of primordial black holes (PBHs) 
is carried out. It is assumed that primordial density fluctuations have power law spectrum and the 
spectral index is constant throughout all scales. It is shown that densities of background radiations 
(y, 7) from PBH evaporations strongly depend on the type of gravitational collapse and on the 
way of taking into account the spread of horizon masses at which PBHs can form. Constraints on 
spectral index values based on PBH evaporation process and on atmospheric and solar neutrino 
■ %■ data are obtained. 
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, I. INITIAL MASS SPECTRA OF PBHS 

The initial mass function of primordial black holes (or, it is better to say, the initial mass spectrum) is determined, 
pL| ' first of all, by the peculiarities of a gravitational collapse near the threshold of a black hole formation. According 
to analytic calculations of seventies a critical size of the density contrast needed for the PBH formation, S c , is 
^vq ' about 1/3. Besides, it was argued that all PBHs have mass roughly equal to the horizon mass at a moment of the 
formation, independently of the perturbation size. However, it was shown recently that near the threshold of black 
hole formation the gravitational collapse behaves as a critical phenomenon In this case the initial mass function 
will be quite different from the analogous function in the analytic calculations Jl],^]. The main feature is that the 
PBH mass formed depends on the size of the fluctuation || , 
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M BH = kM h {8-5 c )~ lk . (1.1) 
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' In this formula Mh is the horizon mass at the time the fluctuation enters the horizon; S c , k, 7 are parameters of 



a concrete model of the critical collapse |3j^]. It is seen from Eq.(Ll) that PBH mass may be arbitrarily small 
independently of the value of Mh- Besides, the value of the critical overdensity, 8 C , in such models is typically ~ 0.7, 
i.e., about a factor of 2 larger than the value found analytically 0. 

The second important ingredient of a PBH initial mass function calculation is a taking into account the spread of 
horizon masses at which PBHs are formed. This problem exists independently of a nature of gravitational collapse. It 
1 is assumed very often that the majority of the black hole formation occurs at the shortest possible scale. In particular, 
the authors of ref. || determined the PBH initial mass function under the assumption that all PBHs form at the same 
horizon mass. The accuracy of such approximation was studied in the work || using the excursion set formalism, and 
it was found that it is good enough, at least in the case of power-law density perturbation spectra with the spectral 
\Ji I index close to 1.2 (i.e., slightly "blue" spectra, with the value of n satisfying observational constraints). 

Recently, in the works of Kim Q and present authors [Q the expression for the PBH mass spectrum (for the 
critical collapse case) was obtained which takes into account the accumulation of PBHs formed at all epochs after 
the beginning of the ordinary radiation dominated era. The derivation is based on the Press- Schechter formalism 
H . Resulting formula for the PBH mass spectrum, taking into account the spread of horizon masses, depends on the 
horizon crossing amplitude er# which is a function of horizon mass Mh- In the present paper we assume a power law 
of primordial density perturbations for all scales, and in this case one has (J1a = 0) 

a H (M h ) = C n S H (M h ), (1.2) 



M eq \ a f M h 



In these formulas 5jj(Af/j) is the non-smoothing horizon crossing amplitude (which is time- independent), M\q is 
the present horizon mass, M eq is the horizon mass at the moment of matter-radiation equality, C n is the numerical 
coefficient arising as a result of the smoothing procedure (the top-hat form of the window-function was used in all 
numerical calculations). The amplitude is normalized on COBE data (see for detailes and references). 
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The expression for the PBH mass spectrum is 

i + 3 [2 
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Here, pi amd Mi are the background density and the horizon mass at ti (ti is the moment of a beginning of growth of 



the density fluctuati ons) . The value of PBH mass Mbh is connected with Mh by the relation (1.1). The upper limit 
of integration in Eq. ( |l.4| ) is determined by the expression 

1/7* 



mm 



M 



BH 



kM, 



(1.5) 



Note, that in ref. the upper limit of integration in the expression for ubh(Mbh) was erroneously put equal to 1. 
All numerical calculations were carried out with the following values of parameters |3j : 



Ik = 0.36 



0.7. 



(1.6) 



The known formula |10(] for the PBH mass spectrum in the standard collapse case can be obtained from Eq.(L4) 
by using the substitutions H 

v 1 / 2 



7 



7 



and the approximate relation 
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Substituting Eqs.(1.7) and (1.8) in Eq.(1.4) one obtains 



tibh{Mbh) = 



W /4 P.M„ 
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In this case the relation between Mbh and Mh is independent on 5: 

M BH =l 1/2 M h . 



(1.7) 
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We carried out in this wor k th e calculations of the neutrino background spectra from PBH evaporations using the 
mass spectrum given by Eq.(1.4) and, in parallel, the mass spectrum from ref. S, 



n B H(M B H) = 
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In this expression the amplitude ub is determined by the same formula_( |l.2| ), as above, but with the constant value 
of Mh, Mh — Mi, in accordance with the assumption of authors of ref. j| that all PBHs form at the smallest scale. 

Our approach has two free parameters: the spectral index n, giving the perturbation amplitude (normalized on 
COBE data on large scales) and ti, the moment of time just after reheating, from which the process of PBHs formation 
started. The value of U is connected with a value of the reheating temperature, 

M p i 



U = 0.301g 



-1/2 



T 2 

RH 
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(<?* ~ 100 is the number of degrees of freedom in the early universe). 

Typical results of PBH mass spectrum calculations are shown on Fig.l. We use, for convenience, the following 
abbreviations : K L, NJ and BK signify, correspondi ngly, the PBH mass spectra calculated by fo rmu lae of Kim and 
Lee 10 1 (Eq. (|l.9| )), Niemeyer and Jedamzik j3| (Eq. Jl.ll )) and Bugaev and Konishchev j?J (Eq.(1.4)). One can see 
from Fig.l that at large values of the spectral index there is a clear difference in a behavior of two mass spectra, NJ 
and BK, based both on a picture of the critical collapse: BK spectrum is less steep at Mbh > Mi. Of course, this 
difference is not so spectacular in cases when the spectral index n is more close to 1 (it is proved by another method 
in ref. §). 

The common feature of NJ and BK spectra is the maximum at Mbh — Mi. Besides, both spectra have same slope 
in the region Mbh < M t . 



2 



II. NEUTRINO BACKGROUND SPECTRA FROM PBHS EVAPORATIONS 



Evolution of a PBH mass spectrum due to the evaporation leads to the approximate expression for this spectrum 
at any moment of time: 

2 

where a accounts for the degrees of freedom of evaporated particles and, strictly speaking, is a function of a running 
value of the PBH mass m. In all our numerical calculations we use the approximation 

a = const = a(M% H x ), (2.2) 

where M^^ x is the value of Mbh m the initial mass spectrum corresponding to the maximum of this spectrum. 
Special study shows that errors connected with such approximation are rather small. 
The expression for the spectrum of the background radiation is H 



= J F(E,z)dlog w (z + l). 



(2.3) 



In this formula f(E, m) is a total instantaneous spectrum of the background radiation (neutrinos or photons) from 
a black hole evaporation. It includes the pure Hawking term and contributions from fragmentations of evaporated 
quarks and from decays of pions and muons (see [^||] for details). The exponential factor in Eq.(^J^) takes into 
account an absorption of the radiation during its propagation in space. The processes of the neutrino absorption are 
considered, in a give n co ntext, in ref. ||. 



In last line of Eq.(2.3) we changed the variable J on z using, for simplicity, the fiat model with = for which 



(2.4) 

n r = (2.4- lO 4 /} 2 )- 1 , ft = 0.67. 



Several examples of z-distributions (the integrands of the integral over z in Eq.( |2.3| )) are shown on Fig. 2. Again 
one can see the rather strong difference of two cases: the existence of a tail of heavy masses in BK spectrum leads 
to a relative enhancement of low z-contributions. The effect becomes more distinct with a rise of the reheating 
temperature. 

The sharp cut-off of all z-distributions near z ~ 10 7 is entirely due to the neutrino absorption ||. The shrinkage of 
z-distributions at large Trh in NJ-case is due to the absence of large masses in the spectrum (PBHs of small masses 
evaporated earlier, and their radiation today is more redshifted). 

Some results of calculations of the neutrino background spectra are shown on Figs. 3-5. The functional form of these 
spectra is qualitatively the same as in the standard collapse case ||: E~ 3 law at large neutrino energies and a flat 
part at low energies with the crossover energy depending on Trh- The only new feature is the sharp steepening of 
the spectra at high values of Trh in NJ case (Fig. 4) connected with the corresponding shrinkage of the z-range. 

Fig. 5 shows the sensitivity of background neutrino intensities to a chosen value of the spectral index n, at a fixed 
Trh- One can see from this figure that the neutrino background flux from PBH evaporations is comparable with the 
atmospheric neutrino flux at ~ 100 MeV if n ~ 1.29 (at T RH ~ 10 9 GeV). Such large values of the spectral index 
are, probably, excluded by the recent large-scale experiments. 

Next figure shows the constraints on the spectral index following from data of neutrino experiments, as a function 
of the reheating temperature. As in our previous work, we used for a description of these constraints the data of the 
Kamiokande atmospheric neutrino experiment jl2| and the experiment on a search of an antineutrino flux from the 
Sun HI (see §] for details). 

For completeness, we show on the same figure the constraints following from extragalactic diffuse gamma ray data; 



part of them (NJ case) is in qualitative agreement with the results of ref. |14|. It is seen that the constraints are much 
more weak in NJ case (at Trh ^ 10 10 GeV). It is again the result of an absence of the tail of large masses in NJ 
spectrum, leading to a shrinkage of the z-distriburions and to a steepening of the background neutrino spectra. 
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III. CONCLUSIONS 



1. The initial concentration of PBH's is larger in the case of the standard picture of the gravitational collapse (as 
compared with the critical collapse case). It is due to relatively small value of the critical overdensity (Sf = 7 = 1/3). 
In the region near the maximum of PBH mass spectra the ratio of intensities is given by the approximate relation 

(KL) si--,? 

^%~^r. (3.i) 

n BH 

Correspondingly, densities of the neutrino and gamma backgrounds produced by PBH's evaporations are smaller and 
spectral index constraints are weaker in the critical collapse case. 

2. The summation over all epochs of PBH's production (i.e., the transition from NJ to BK case) leads to an 
appearing of the tail of large PBH masses in the PBH mass spectrum. This tail is especially sufficient if the spectral 
index n is relatively large. 

3. In scenarios with high reheating temperatures, Trh ~ 10 9 — 10 10 GeV, the summation over all epochs leads to a 
strong enhancement of a high energy part of the background neutrino spectra. 

4. Spectral index constraints following from the comparison of neutrino background predictions with existing data 
of neutrino experiments are rather weak if the purely power law of primordial density fluctuations is assumed for all 
scales (and, correspondingly, the normalization on COBE data is used). However, it does not mean that neutrino 
background intensities from PBH's evaporations cannot be noticeable. If, for example, the spectrum of primordial 
density fluctuations is a combination of two simple power law spectra |6j, i.e., 

<J H (M h ) ~ , M h < M hc < M eq ; 

(3.2) 

(T H (M h ) ~ , M eq > M h > M hc , 

then the small value of rii ~ 1 (n; — 1 ~ 0) which follows, in particular, from COBE data, does not contradict with the 
possibility that n s — 1 is large and, correspondingly, an at small scales is also large. If Mh c and n; are known from 
somewhere, one can easily obtain from the neutrino data the constraints on the small scale spectral index n s , in the 
same manner as it is done above for the all scale spectral index n. Several cosmological models of such kind, which 
predict a large density fluctuation amplitude just on small scales (while the amplitude on large scales is constrained 
by the small CMBR anisotropy) appeared recently (see, e.g., |^5| , |l6| ). 
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FIG. 1. Examples of PBH mass spectra for two values of Trh- Solid lines: KL spectra [10], dotted lines: NJ spectra [3], 
dashed lines: BK spectra [7]; n = 1.30 in all cases. 
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FIG. 2. Red-shift distributions (integrands of the expression (2.3) for the background spectrum) for neutrino energy 
E — 1CP 1 GeV (solid lines) and E = 10~ 2 GeV (dashed lines), for PBH mass spectrum (left column) and BK mass spec- 
trum (right column), for three values of Trh', n = 1.30 in all cases. The dimension of F(z) is (s~ 1 cm~ 2 sr~ 1 GeV~ 1 ). 
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FIG. 3. Neutrino background spectra for three values of Trh, for BK mass spectrum; n = 1.30. 
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FIG. 5. Neutrino background spectra for three values of the parameter n, for BK mass spectrum. Dashed line represents the 
theoretical atmospheric neutrino spectrum for Kamiokande site (averaged over all directions) [11]. 
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FIG. 6. Constraints on the spectral index n as a function of the reheating temperature Trh- Solid and dotted lines: BK 
and NJ mass spectra, atmospheric u e and solar v e experiments. Dashed and dot-dashed lines: BK and NJ mass spectra, 
extragalactic diffuse gamma-ray background data. 
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